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Abstract. We study the procedure of the reconstruction of phantom-scalar field 

potentials in two-field cosmological models. It is shown that while in the one-field case 

C^ ' the chosen cosmological evolution defines uniquely the form of the scalar potential, in 

the two-field case one has an infinite number of possibilities. The classification of a large 
class of possible potentials is presented and the dependence of cosmological dynamics 
on the choice of initial conditions is investigated qualitatively and numerically for two 
particular models. 
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1. Introduction 

The discovery of cosmic acceleration \\\ has stimulated the construction of a class of 
dark energy models [2],[3l[4j describing this effect. Notice that the cosmological models 
based on scalar fields were considered long before the observational discovery of cosmic 
acceleration [5]. The dark energy should possess a negative pressure such that the 
relation between pressure and energy density w is less than —1/3. The present day 
value of the parameter w is close to —1 (cosmological constant) but some observations 
indicate that the value of this parameter slightly less than —1 provides the best fit. 
The corresponding dark energy has been named phantom dark energy [6j. According to 
some authors, the analysis of observations permits the existence of the moment when 
the universe changes the value of the parameter w from w > —1 to w < —1 [TllH]- 
This transition is called "the crossing of the phantom divide line". The most recent 
investigations have shown that the phantom divide line crossing is still not excluded by 
the data [9]. 

It is easy to see, that the standard minimally coupled scalar field cannot give rise 
to the phantom dark energy, because in this model the absolute value of energy density 
is always greater than that of pressure, i.e. |ti;| < 1. A possible way out of this situation 
is the consideration of the scalar field models with the negative kinetic term. Thus, 
the important problem arising in connection with the phantom energy is the crossing of 
the phantom divide line. The general belief is that while this crossing is not admissible 
in simple minimally coupled models its explanation requires more complicated models 
such as "multifield" ones or models with non-minimal coupling between scalar field and 
gravity (see e.g. pOj). 

In preceding papers some of us [HI [12] described the phenomenon of the change 
of sign of the kinetic term of the scalar field implied by the Einstein equations. It was 
shown that such a change is possible only when the potential of the scalar field possesses 
some cusps and, moreover, for some very special initial conditions on the time derivatives 
and values of the considered scalar field approaching to the phantom divide line. At the 
same time, two-field models including one standard scalar field and a phantom field can 
describe the phenomenon of the (de)-phantomization under very general conditions and 
using rather simple potentials ^3\, [H]. In the present paper we would like to attract 
attention to a drastic difference between two- and one-field models. 

The reconstruction procedure of the scalar field potential models is well-known 
[ISlIinillTlIIHlIiniEOlEIlES]. We recapitulate it briefly. The cosmological evolution in 
a flat Friedmann model with the metric 

ds^ = dt'-a^{t)df (1) 

can be given by the time evolution of the Hubble parameter h{t) = ^, which satisfies 
the Friedmann equation 

h^ = e, (2) 
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where e is the energy density and a convenient normalization of the Newton constant is 
chosen. Differentiating equation ^ and using the energy conservation equation 



e 



-3hie + p), (3) 



where p is the pressure, one comes to 

h = --ie + p). (4) 

If the matter is represented by a spatially homogeneous minimally coupled scalar field, 
then the energy density and the pressure are given by the formulae 

s=l<P' + V{<f>), (5) 

P=^0'-n0), (6) 

where V{(j)) is a scalar field potential. Combining equations ([2]), (l4|), ([5]), (JHI we have 

y=^ + h^ (7) 

and 

02 = -h. (8) 

Equation ((71) represents the potential as a function of time t. Integrating equation 
(IH|) one can find the scalar field as a function of time. Inverting this dependence we 
can obtain the time parameter as a function of (j) and substituting the corresponding 
formula into equation ([7]) one arrives to the uniquely reconstructed potential V"(0). It is 
necessary to stress that this potential reproduces a given cosmological evolution only for 
some special choice of initial conditions on the scalar field and its time derivative [iTlfTI] . 
Below we shall show that in the case of two scalar fields one has an enormous freedom 
in the choice of the two-field potential providing the same cosmological evolution. This 
freedom is connected with the fact that the kinetic term has now two contributions. In 
order to examine the problem of phantom divide line crossing we shall be interested here 
in the case of one standard scalar field and one phantom field ^ , whose kinetic term 
has a negative sign. In this case the total energy density and pressure will be given by 

e=\<P'-li' + V{<P,0, (9) 

p=^0'-^^'-n0,O- (10) 

The relation ((71) expressing the potential as a function of t does not change in form, but 
instead of equation (IHj) we have 

^'-e = -lh. (11) 

Now, one has rather a wide freedom in the choice of the time dependence of one of two 
fields. After that the time dependence of the second field can be found from equation 
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f lTTl) . However, the freedom is not yet exhausted. Indeed, having two representations 
for the time parameter t as a function of or ^, one can construct an infinite number of 
potentials V^(0, using the formula (jTl) and some rather loose consistency conditions. 
It is rather difficult to characterize all the family of possible two-field potentials, 
reproducing given cosmological evolution h{t). In the present paper, we describe some 
general principles of construction of such potentials, then we consider some concrete 
examples. 

2. The system of equations for two-field cosmological model 

The system of equations, which we study contains ([7]) and ( fTTl) and two Klein-Gordon 
equations 

^ + 3/.^ +^5^ = 0. (12) 

{■+3^-^ = 0. (13) 

From equations flT2|l and (IT3l) we can find the partial derivatives — -^^ and — 1|^ as 
functions of time t. The consistency relation 

is respected. 

Before starting the construction of potentials for particular cosmological evolutions, 
it it useful to consider some mathematical aspects of the problem of reconstruction of a 
function of two variables in general terms. 

2.1. Reconstruction of the function of two variables, which in turn depends on a third 
parameter 

Let us consider the function of two variables -F(x, y) defined on a curve, parameterized 
by t. Suppose that we know the function F{t) and its partial derivatives as functions of 

t: 

F{xit),y{t)) = F{t), (15) 

dFixit),yit)) dF 

" l^yt), (16) 



dx dx 

dF{x{t),y{t)) OF 



it). (17) 



dy dy 

These three functions should satisfy the consistency relation 

OF OF 

F(t) = — (t)x + — (t)y. (18) 
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As a simple example we can consider the curve 

x(t)=t, y{t) = t\ (19) 
while 

Fit) = t\ (20) 

dF , , 

1% = '• '''' 

and equation ( ITSll is satisfied. 

Thus, we would like to reconstruct the function F{x,y) having explicit expressions 
in right-hand side of equations fITSil- ( flTll . This reconstruction is not unique. We shall 
begin the reconstruction process taking such simple ansatzes as 

F,{x,y) = G^ix)+H,iy), (23) 

F2{x,y) = G2ix)H,{y), (24) 

F3{x,y) = {G3{x) + Hs{y)r. (25) 



(26) 
(27) 



The assumption (12311 immediately implies 
dFi _ dGi 

dx dx 
dFi _ dHi 
dy dy 

Therefrom one obtains 

-^{t{x'))dx\ (28) 

H,{y) = l'^{t{y'))dy'. (29) 

Hence 

^{t{x'))dx' + J -^{t{y'))dy'. (30) 

For an example given by equations (fT9l)- (|22l) the function Fi{x,y) is 

/x py px py 

t'^{x')dx' + / t{y')dy' = / x'^dx' ± /' y^dy. (31) 

Explicitly 

\ {x' + 2y^'^) , x>0, 
Fiix,y)= { f (32) 

- {x' - 2|/3/2) , x<0. 
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Similar reasonings give for the assumptions (l24l) .and (i25ll correspondingly 

F^ix, y) = exp 11 (}^) {t{x))dx + j (-^^^ {t{y))dy^ , (33) 

For our simple example (IT9ll - (l22l) the functions F2{x,y) and F^{x,y) have the form 

F2{x,y)=xy, (35) 

'1 



F?.{x,y) 



{x^l^ + 2I/3/2-) 



i (x3/"(-)3/"2y3/2a) 



x>0, 

(36) 
x<0. 



Thus, we have seen that the same input of "time" functions fl20l)-fl22ll on the curve (ITOD 
produces quite different functions of variables x and y. 

Naturally, one can introduce many other assumptions for reconstruction of -F(x, y). 
For example, one can consider linear combinations of x and y as functions of the 
parameter t and decompose the presumed function F as a sum or a product of the 
functions of these new variables. 

Now we present a way of constructing the whole family of solutions starting from a 
given one. Let us suppose that we have a function Fo{x,y) satisfying all the necessary 
conditions. Let us take an arbitrary function 

/(.,.)^/(||). (37) 

which depends only on the ratio t{x)/t{y). We require also 

/(1) = 1, /(1) = 0, (38) 

i.e. the function reduces to unity and its derivative vanishes on the curve {x(t),y{t)). 
Then it is obvious that the function 

F{x,y) = Foix,y)f(^-^) (39) 

is also a solution. This permits us to generate a whole family of solutions, depending 
on a choice of the function /. Moreover, one can construct other solutions, adding to 
the function F{x,y) a term proportional to (t{x) — t{y)Y. 

2.2. Cosmological applications, an evolution "Bang to Rip" 

To show how this procedure works in cosmology, we consider a relatively simple 
cosmological evolution, which nevertheless is of particular interest, because it describes 
the phantom divide line crossing. Let us suppose that the Hubble variable for this 
evolution behaves as 

Ht) = -, r, (40) 

^^ t{tR-t) ^ ' 
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where A is a positive constant. At the beginning of the cosmological evolution, when 
t — > the universe is born from the standard Big Bang type cosmological singularity, 
because h{t) ~ 1/t. Then, when t — > t^, the universe is superaccelerated, approaching 
the Big Rip singularity h{t) ~ l/{tR — t). Substituting the function (l40ll and its time 
derivative into equations. (fTTl) and (jTl) we come to 

For convenience we choose also the parameter A as 

^ = I , (43) 

Then, 

and 

Let us consider now a special choice of functions (f){t) and ^(t) used alreadjljin [T3] . 

0(t) = — arctanhW , (46) 

3 \ tR 

The derivatives of the potential with respect to the fields (j) and ^ could be found from 
the Klein-Gordon equations ( fT2ll and (IT3ll : 



We can obtain also the time parameter as a function of or ^: 

^(^) - cosh-(-30/4) ' ^^°^ 

^^^^ " exp(-3e/v^) + l' ^^^^ 

Now we can make a hypothesis about the structure of the potential V{^, (j)): 

V2{iA) = G{m{<P)- (52) 

X Notice that the origin of two scalar fields has been associated in [TSj with a non-Hermitian complex 
scalar field theory and there a classical solution was found as a saddle point in "double" complexification. 
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Applying the technique described in the subsection A, we can get G{^): 

and 

G(0=exp(3^0- (54) 

To find H{(j)) one can use the analogous direct integration, but we prefer to implement 
a formula 

-- ^W*> (B5) 



(B6) 





"'^'-G((m)y 


which gives 






Him) = If, 


and hence, 








Finally, 






V2{^A) = ^cosh6(-3</./4)exp(3V20• 


Here we have reproduced the potential studied in |1.Sj. 


Making 


the choice 




^i(e,0) = G(O + i^(0), 


we derive 





(57) 



(58) 



(59) 



^.(e,0)-^ 



_le"3?/v^ + 3v/2e + 2e3^/^ + -e^^/^ 



sinh^(-30/4) + sinh^(-30/4) - 1 , , ., , ,,/ 
^ — ^' \^ ^^ \^ ^' ' + lnsinh^(-30/4) 



sinh^(-30/4) ' v -^/ -. • K ) 

Now, we can make another choice of the field functions 0(t) and ^(t), satisfying the 
condition f l4Tl) : 

0(t) = ^ln^, (61) 

i{t) = ^arctanh J^. (62) 

The time parameter t is a function of fields is 

t(0) = ^^^ , (63) 

^^ exp(-30/v^) + 1 ^ ^ 
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tR 



m 



tanh2(3e/4)' 



9 

(64) 



Looking for the potential as a sum of functions of two fields as in equation ([59 
after lengthy but straightforward calculations we come to the following potential: 



Vi{^, 



3tl 



1 + - sinh^(3^/4) + 3sinh2(3^/4) 



3sinh^(3^/4) 



+ 21nsinh2(3^/4) + - exp(- 30/^2) - 3^20 
o 

-2exp(30/V2) - -exp(0/V2) 



(65) 



Similarly for the potential designed as a product of functions of two fields (l52!l we 
obtain 



V2{^, (t)) = TrT sinh2(3^/4) cosh^(3^/4) exp(-3V20). 



9t 



(66) 



R 



One can make also other choices of functions (pit) and ^(t) generating other 
potentials, but we shall not do it here (see [23]), concentrating instead on the qualitative 
and numerical analysis of two toy cosmological models described by potentials (l65l) and 
(El. 



3. Analysis of cosmological models 



It is well known [24] that for the qualitative analysis of the system of cosmological 
equations it is convenient to present it as a dynamical system, i.e. a system of first- 
order diff'erential equations. Introducing the new variables x and y we can write 





' = X, 






^ =y, 




< 


X = -3sign(/i)xA/ Y 


-^ + ne,0) 




^ y = -3sign(/i)yy ^ 


-f + ne,0) 


ice that the refiection 




x-^-x, 




y^ -y, 




t 


— > -t 





dV 
dV 

ar 



(67) 



(68) 

transforms the system into one describing the cosmological evolution with the opposite 
sign of the Hubble parameter. The stationary points of the system ( l67l l are given by 



dV dV 

x = 0, , = 0, — = 0, — 



0. 



(69) 
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In this subsection we shall analyze the cosmological model with two fields - standard 
scalar and phantom, described by the potential (l65l ). For this potential the system of 
equations (l67l) reads 



X, 



^ 

X 



— 3sign(/i)x 



- ^ + 



9ti 



■sinh2(3^/4)cosh^(3^/4)e^ 



^ZV2(t> 



+ ^ sinh2(3^/4) cosh^(3e/4)e- 



3V2(/> 



(70) 



ij = -3sign{h)yJ'^ - T + ^ sinh2(3^/4) cosh^(3^/4)e-3v^<^ 



sinh(3g/4)cosh'^(3g/4) ri 



3ti 



[| cosh2(3e/4) + sinh2(3^/4)] e-'^^^. 



It is easy to see that there are stationary points 

= 00, e = 0, X = 0, y = 0, (71) 

where 0o is arbitrary. For these points the potential and hence the Hubble variable 
vanish. Thus, we have a static cosmological solution. We should study the behavior of 
our system in the neighborhood of the point (FTTl ) in linear approximation: 



X 

y 



y, 



I _C p-3v/2(/>o 



(72) 



One sees that the dynamics of in this approximation is frozen and hence we 
can focus on the study of the dynamics of the variables ^,y. The eigenvalues of the 
corresponding subsystem of two equations are 

g-30o/v^ 

Ai,2 = T^^T • (73) 

These eigenvalues are real and have opposite signs, so one has a saddle point in the 
plane (^, y) and this means that the points (iTTll are unstable. 

One can make another qualitative observation. Freezing the dynamics of ^ 
independently of 0, namely choosing y = 0,^ = 0, which implies also y = C, = 0, 
one has the following equation of motion for 0: 

+ 3/i0 = 0. (74) 

Equation (rf4l) is nothing but the Klein-Gordon equation for a massless scalar field on 
the Friedmann background, whose solution is 

V2. 



m 



Int, 



and which gives a Hubble variable 

1 
3t' 



h{t) 



(75) 



(76) 
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This is an evolution of the flat Friedmann universe, filled with stiff matter with the 
equation of state p = e. It describes a universe, born from the Big Bang singularity and 
infinitely expanding. Naturally, for the opposite sign of the Hubble parameter, one has 
the contracting universe ending in the Big Crunch cosmological singularity. 




Figure 1. An example of section of the 4D phase space obtained with numerical 
calculations. We can see a series of trajectories corresponding to different choices of 
the initial conditions. Every initial condition in the graphic is chosen on the dashed 
"ellipse" centered in the origin and is emphasized by a colored dot. A "saddle point-like" 
structure of the set of the trajectories clearly emerges. 



Now, we describe some results of numerical calculations to have an idea about 
the structure of the set of possible cosmological evolutions coexisting in the model 
under consideration. We have carried out two kinds of simulations. First, we have 
considered neighborhood of the plane y, ^ with the initial conditions on the field such 
that 0(0) = 0,0(0) = (see Figure [T|). The initial conditions for the phantom field 
were chosen in such a way that the sum of absolute values of the kinetic and potential 
energies were fixed. Then, running the time back and forward we have seen that the 
absolute majority of the cosmological evolutions began at the singularity of the "anti-Big 
Rip" type (Figure [2]). Namely, the initial cosmological singularity were characterized 
by an infinite value of the cosmological radius and an infinite negative value of its 
time derivative (and also of the Hubble variable). Then the universe squeezes, being 
dominated by the phantom scalar field ^. At some moment the universe passes the 
phantom divide line and the universe continues squeezing but with h < 0. Then it 
achieves the minimal value of the cosmological radius and an expansion begins. At 
some moment the universe undergoes the second phantom divide line crossing and 
its expansion becomes super-accelerated culminating in an encounter with a Big Rip 
singularity. Apparently this scenario is very different from the standard cosmological 
scenario and with its phantom version Bang-to-Rip, which has played a role of an 
input in the construction of our potentials. The second procedure, which we have used 
is the consideration of trajectories close to our initial trajectory of the Bang-to-Rip 
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Figure 2. Typical behavior of the Hubble parameter for the model I trajectories. 
The presence of two stationary points (namely a maximum followed by a minimum) 
indicates the double crossing of the phantom divide line. Both the initial and final 
singularities are characterized by a Big Rip behavior, the first is a contraction and the 
second is an expansion. 



type. The numerical analysis shows that this trajectory is unstable and the neighboring 
trajectories again have anti-Big Rip - double crossing of the phantom line - Big Rip 
behavior described above. However, it is necessary to emphasize that a small subset of 
the trajectories of the Bang-to- Rip type exist, being not in the vicinity of our initial 
trajectory. 



3.2. Model II 

In this subsection we shall study the cosmological model with the potential f l58ll . Now 
the system of equations (l67ll looks like 



x, 



X 



+ 2^cosh6(30/4)e3v^«. 
Notice that the potential (l58ll has an additional reflection symmetry 



-3sign(/i)xy/f - \ + ^cosh^(30/4) exp{3v^e} 
-7^ cosh^(30/4) sinh(30/4)e3v^«, 

'-R 

-3sign(/i)|/./f - f + q^ cosh^(30/4) exp{3v^0 



(77) 



(78) 
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This provides the symmetry with respect to the origin in the plane (0, x). The system 
(jTTll has no stationary points. However, there is an interesting point 

= 0, X = 0, (79) 

which freezes the dynamics of and hence, permits to consider independently the 
dynamics of ^ and y, described by the subsystem 

y= -3sign(/i)yy/-^ + ^e3v^€ + f^e^v^?. ^^^> 

Apparently, the evolution of the universe is driven now by the phantom field and is 
subject to superacceleration. 

In this case the qualitative analysis of the differential equations for ^ and y, 
confirmed by the numerical simulations gives a predictable result: being determined 
by the only phantom scalar field the cosmological evolution is characterized by the 
growing positive value of h. Namely, the universe begins its evolution from the anti- 
Big Rip singularity [h = — oo) then h is growing passing at some moment of time the 
value h = (the point of minimal contraction of the cosmological radius a(t)) and then 
expands ending its evolution in the Big Rip cosmological singularity {h = +oo). 

Another numerical simulation can be done by fixing initial conditions for the phan- 
tom field as ^(0) = 0,y{0) = (see Figure [3]). Choosing various values of the initial 
conditions for the scalar field 0(O),a;(O) around the point of freezing = 0,a; = we 
found two types of cosmological trajectories: 

1. The trajectories starting from the anti-Big Rip singularity and ending in the Big Rip 
after the double crossing of the phantom divide line. These trajectories are similar to 
those discussed in the preceding subsection for the model I. 

2. The evolutions of the type Bang-to-Rip. 

Then we have carried out the numerical simulations of cosmological evolutions, 
choosing the initial conditions around the point of the phantomization point with the 
coordinates 

0(0) = 0, 

/ X 4 1 

4(0) = -arctanh^p, 
3 v2 

1/(0) = ^. (81) 

This analysis shows that in contrast to the model I, here the standard 
phantomization trajectory is stable and the trajectories of the type Bang-to-Rip are 
not exceptional, though less probable then those of the type anti-Big Rip to Big Rip. 
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Figure 3. An example of section of the 4D phase space for the model II obtained with 
numerical calculations. We can see a family of different trajectories corresponding 
to different choices of the initial conditions. Every initial condition in the graphic is 
chosen on the dashed "ellipse". The origin represents the point of freezing. 



4. Conclusion 



We have considered the problem of reconstruction of the potential in a theory with 
two scalar fields (one standard and one phantom) starting with a given cosmological 
evolution. It is known ( see e.g. (IH [HI [iTj) that in the case of the only scalar field 
this potential is determined uniquely as well as the initial conditions for the scalar 
field, reproducing the given cosmological evolution. Changing the initial conditions, 
one can find a variety of cosmological evolutions, sometimes qualitatively different from 
the "input" one (see e.g. tllj)- In the case of two fields the procedure of reconstruction 
becomes much more involved. As we have shown here, there is a huge variety of different 
potentials reproducing the given cosmological evolution (a very simple one in the case, 
which we have explicitly studied here). Every potential entails different cosmological 
evolutions, depending on the initial conditions. 

It is interesting that the existence of different dynamics of scalar fields corresponding 
to the same evolution of the Hubble parameter h{t) can imply some observable 
consequences connected with the possible interactions of the scalar fields with other 
fields. Indeed, in this case, the time dependence of the scalar fields considered above 
can directly affect physically observable quantities. We are going to consider this topic 
elsewhere 1251. 
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